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• Abstract. I describe a model of energetic gamma ray and neutrino 

emission in blazars which is consistent with most of the gamma-ray ob- 
servations, and use this model to estimate the diffuse intensity of high 
energy neutrinos from blazars. 
*Accretion Phenomena and Related Outflows, IAU Colloq. 163, 1996 
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Introduction 



. The second EGRET catalog of high-energy gamma-ray sources (Thompson 1995) 

\^ | contains 40 high confidence identifications of AGN, and all appear to be blazars. 

Clearly, the gamma ray emission is associated with the jet. Several of these 
C""-- AGN show gamma-ray variability with time scales of ~ 1 day (Kniffen 1993), 

and this places an important constraint on the models. 
\ Most previous work on gamma ray emission in AGN jets involved elec- 

tron acceleration and inverse Compton scattering. The proton blazar model 
was proposed originally by Mannheim and Biermann (1992) and developed by 
Mannheim (1996). In this model both protons and electrons are accelerated, 
and protons interact via pion photoproduction and Bethe-Heitler pair produc- 
£^ ■ tion with synchrotron radiation produced by the electrons. Here I consider a 

\ model in which protons are accelerated and interact directly with accretion disk 

radiation. 
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■ 2. The model 

a i 

If the protons interacted with matter to produce the gamma rays on a timescale 
consistent with the observed variability, then the proton number density in the 
blob must be nn > 10 9 cm -3 , assuming a Doppler factor of ~ 20. Such a 
high density would not normally be expected in AGN jets. Interactions with the 
radiation field are therefore necessary if protons are accelerated, but a sufficiently 
high density of target photons for proton interactions would present problems 
for gamma ray escape if the radiation were isotropic. In the present work, 
the gamma-ray emission region (referred to as "the blob") is placed well above 
the accretion disk, and any scattering by clouds or coronae is assumed to be 
negligible, so that the radiation impinging on the blob is highly anisotropic. 
I use a spectrum which is a modification of the standard thin accretion disk 
spectrum (Shakura and Sunyaev 1973) in which the maximum temperature is 
lower than usual to account for the absence of the inner part of the disk due 
to the presence of the jet. This is only an approximation, and does not include 
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the soft X-ray emission. Calculations are in progress using a more appropriate 
spectrum and will be reported elsewhere. However, the results for GeV gamma- 
rays and high energy neutrinos are not expected to differ by very much from 
those presented here. 

From variability arguments, the size of the blob can not be larger than 
cAt' ygx (primed quantities refer to the jet frame). This means that the region 
in which the protons are accelerated and trapped must be smaller than 0.02 pc, 
and proton gyroradii, r' g , at the highest energies must be less than 0.01 pc. The 
maximum acceleration rate for any acceleration mechanism is (dE'/dt') max = 
ec 2 B' where B' is the magnetic field in the blob, but for realistic acceleration 
mechanisms the acceleration rate will be at least a factor of 10 lower. During 
acceleration, protons will also suffer synchrotron losses and we require that these 
be less than the energy losses for pion photoproduction at all energies up to 
the maximum energy (which is determined by pion photoproduction). For a 
magnetic field of 10 gauss all these conditions are met, giving a maximum jet- 
frame energy of E' max ~ I- 5 x 1° 10 GeV. 

Protons may escape from the blob before interacting, and I assume this is 
simply due to diffusion to the boundary of the blob where they freely escape. For 
a reasonable diffusion coefficient, most protons with energies below the threshold 
for pion photoproduction will escape without interacting. 

3. Interaction of protons and the pair-synchrotron cascade 

We inject an E'~ 2 spectrum of protons, continuously randomize their directions 
in the jet frame, and allow them to interact with radiation. The interactions 
are treated as described in Szabo and Protheroe (1994) and all pions, and their 
decay products, are assumed to be produced travelling initially in the direction 
of the interacting proton. 

The gamma rays from 7r° decay may interact with photons from the ac- 
cretion disk while inside the blob, or outside the blob, or may escape from the 
AGN. The probability of escaping from the blob and escaping from the AGN are 
calculated for various angles to the jet axis. The gamma rays that pair produce 
inside the blob will initiate a pair-synchrotron cascade. The electrons produced 
inside the blob by photon-photon pair production of gamma-rays from 7r° de- 
cay, and electrons from from tt^ decay, will synchrotron radiate in the ambient 
magnetic field and produce the first generation of synchrotron photons in the 
cascade. I then treat the photon-photon pair production interactions of this first 
generation of synchrotron gamma-rays in the same way as gamma rays from 7r° 
decay, the resulting electron positron pairs producing the second generation of 
synchrotron photons in the cascade. This procedure is repeated for several gen- 
erations. The resulting gamma ray and neutrino fluxes emerging from the AGN 
at angles between 2° and 4° with respect to the jet axis are shown in Fig. 1. 

4. Neutrino background from unresolved blazars 

I use the luminosity function of gamma-ray loud AGN (Chiang 1995) and inte- 
grate over redshift and luminosity to estimate the diffuse background shown in 
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Figure 1. Average gamma ray and neutrino spectra for viewing an- 
gles 2° - 4° with respect to the jet axis resulting from jry interactions 
and subsequent cascading. Full hisogram shows the total neutrino out- 
put (i/fj, + v lJL + Ve + v e ), dotted hisogram shows the 7r° gamma spectrum 
on production, and dashed, dot-dashed, and dot-dot-dot-dashed curves 
show 1st, 2nd and 3rd generation synchrotron spectra on production. 
Full curve shows the emerging gamma ray spectrum. 



Fig. 2. Also shown is the result of Mannheim (1995) which, given the differences 
in the models, is in surprisingly good agreement with the present work. The 
reason for this is that despite the model differences, the most important param- 
eters have similar values (e.g. we inject the same proton spectrum, have similar 
maximum energies, similar magnetic fields, and similar maximum soft photon 
energies) . 

The predicted neutrino intensities may be detectable with the generation 
of neutrino telescopes currently being constructed (see Gaisser et al. 1995 for a 
review). For the present work, the number of muon events per year in a detector 
with an effective collecting area of 2 x 10 4 m 2 would be ~ 10 above 1 TeV, and 
~ 6 above 10 TeV, compared with the atmospheric neutrino background of 
~ 225 and 5 per year respectively (Hill 1996). 

As a final note, I mention that one of the EGRET blazars, Markarian 421, 
which was previously observed at TeV energies (Punch et al., 1992), very recently 
had two violent outbursts at TeV energies with a variability timescale of ~ 1 h 
(Gaidos et al. 1996). Such rapid variability has yet to be adequately explained 
by any of the models, including that described in the present work. 
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Figure 2. The expected diffuse + intensity compared with the 
atmospheric neutrino intensity (Lipari 1993) (vertical hatched band: 
the upper curve corresponds to zenith angle 9 = 90° and the lower 
curve corresponds to 9 = 0°). Also shown are possible minimum and 
maximum backgrounds due to charm production (dashed lines). Dot- 
ted curve shows the result of Mannheim (1995). 

References 

Chiang, J., et al., ApJ, 452 156 (1995). 
Gaidos, J. A., et al., Nature in press (1996). 

Gaisser, T.K., Halzen, F., and Stanev, T., Physics Reports, 258, 173 (1995). 

Hill, G., Astroparticle Phys. in press (1996). 

Kniffen, D.A., et al., ApJ, 411 133 (1993). 

P. Lipari, Astroparticle Phys. 1 195 (1993). 

Mannheim, M., and Biermann, P.L., A&A221, 211 (1989). 

Mannheim, K., Astroparticle Phys. 3 295 (1995). 

Punch, M., et al. Nature 358, 477 (1992). 

Shakura, N.I., and Sunyaev, R.A., A&A, 24 337 (1973). 

Szabo A. P. and Protheroe R.J. Astroparticle Phys. 2 375 (1994). 

Thompson, D.J., et al., ApJS, 101 259 (1995). 



4 



